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CN ' Abstract 

m . 

O . We extend to all orders in perturbation theory a method to calculate 



supersymmetry-breaking effects by analytic continuation of the renormal- 

ization group into superspace. A central observation is that the renormal- 

Ph! ized gauge coupling can be extended to a real vector superfield, thereby 

O-i. including soft breaking effects in the gauge sector. We explain the rela- 

D ■ 

r-| . tion between this vector superfield coupling and the "holomorphic" gauge 

coupling, which is a chiral superfield running only at 1 loop. We consider 
^ ! these issues for a number of regulators, including dimensional reduction. 

^ i With this method, the renormalization group equations for soft supersym- 

metry breaking terms are directly related to supersymmetric beta functions 
and anomalous dimensions to all orders in perturbation theory. However, 
the real power of the formalism lies in computing finite soft breaking ef- 
fects corresponding to high-loop component calculations. We prove that 
the gaugino mass in gauge-mediated supersymmetry breaking is "screened" 
from strong interactions in the messenger sector. We present the complete 
next-to-leading calculation of gaugino masses (2 loops) and sfermion masses 
(3 loops) in minimal gauge mediation, and several other calculations of phe- 
nomeno logical relevance. 
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1 Introduction 



Recently there has been a great deal of interest in building models in which supersym- 
metry breaking is communicated to the observable particles through renormalizable 
interactions A common feature of these models is that supersymmetry breaking 
occurs in the masses of "messenger" fields in the form 

M = MsusY + 5M, (1.1) 

where Msusy is a supersymmetric mass term, and 5M breaks supersjTiimetry. In 
most models of this kind constructed to date 5M <C Msusy, and so the messenger 
threshold is approximately supersymmetric. Integrating out the messenger fields gives 
rise to supersymmetry breaking in the low-energy effective lagrangian below the scale 
M. A large amount of work has already been done on the calculation of the su- 
persymmetry breaking effects from various types of interactions 0, 0]. In Ref. it 
was shown how to compute the leading low-energy supersymmetry breaking effects 
in a large class of models using only one-loop renormalization group (RG) equations 
and tree-level matching, while direct calculations of the same quantities require the 
evaluation of 1- and 2-loop graphs. 

The starting point of Ref. is the observation that since the messenger threshold 
is approximately supersymmetric, one can use a formalism where all couplings and 
masses are treated as superfields, and the SUSY breaking terms correspond to non- 
zero ^-dependent spurion components of the couplings. In this framework, it is not 
hard to see that leading-log effects that are determined by the RG in the SUSY limit 
are related to finite SUSY-breaking effects. For example, the RG can be used to 
compute corrections of the form (In M)/(167r^), where M is a threshold mass. If M 
is a superfield, then this contribution has a SUSY-breaking component 

J_ InML,,, = (1.2) 
167r2 167r2 M\^ ' ^ ' 

which contains a loop factor, but no logarithm. Effects of this type therefore corre- 
spond to finite loop effects that are not related to an RG calculation in components. 

A simple power-counting argument can be used to show that in gauge-mediated 
models the leading SUSY-breaking terms in the low-energy effective lagrangian arise 
from this sort of threshold dependence in the dimensionless couplings. This allows one 
to compute 1- and 2- loop SUSY breaking effects using the 1-loop RG equations and 
tree-level matching, analytically continued into superspace. In Ref. this technique 
was used to reproduce known results in a much simpler way, and also to derive new 
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phenomenologically interesting results that would be much more difficult to compute 
directly. 

In this paper, we extend the analysis of Ref. to higher orders in perturbation 
theory. One motivation for this is to define an unambiguous procedure to perform 
the analytic continuation into superfield beyond one loop. We show that the gauge 
coupling is naturally extended to a real superfield that is not the sum of a chiral and an 
antichiral superfield. The O'^O'^ component of the real gauge superfield plays a crucial 
role in reproducing the correct behavior of perturbation theory. Another motivation 
for this is to obtain new results of interest for testing models in the literature. In 
particular, we are able to compute gaugino, squark, and slepton masses in gauge- 
mediated models at the next-to-leading order in perturbation theory. Our result 
corresponds to an explicit calculation of 2- and 3-loop Feynman diagrams. One of 
our results is that the gaugino masses in gauge-mediated models are "screened" from 
corrections from the SUSY-breaking sector up to 4 loops. This implies that the gauge- 
mediation relations are preserved up to corrections of order g^y^/ {IQn'^Y ~ 10"'' even 
if the SUSY-breaking (or messenger) sector is strongly coupled. We also compute 
other interesting effects, like the gaugino masses in "mediator" models @|, the gauge- 
mediated effective potential induced along classically flat directions, both for D fiat 
directions (2 loops) as well as for the scalar partner of the axion (3 loops). 

This paper is organized as follows. In Section 2, we give a definition of renor- 
malized coupling constants that can be viewed as superfield spurious to all orders in 
perturbation theory. We use as examples specific theories that allow simple super- 
symmetric regulators. In Section 3, we discuss this prescription in the case in which 
the theory is regulated using dimensional reduction. We also show that extending the 
couplings to superfields automatically selects the so-called DR' scheme for the soft 
terms. In Section 4, we use our technique to prove the gaugino screening result men- 
tioned above, and compute gluino, squark, and slepton masses in gauge mediation at 
the NLO. We also extend our results to D-term breaking of SUSY, and derive the 
gaugino mass in "mediator" models. In Section 5, we compute some other interesting 
SUSY-breaking effects in gauge-mediated theories. Section 6 summarizes our main 
results and contains our conclusions. 

2 Renormalized Coupling Constants as Superfields 

The main tool of our approach is the use of renormalization schemes in which the 
renormalized coupling constants can be treated as superfields. Much of our discussion 
can be viewed as a restatement of the insights of Shifman and Vainshtein ||^ in 
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the framework of renormalized perturbation theory. However, we will generalize the 
method to include supersymmetry breaking effects. For gaugino masses and A terms, 
this was first done in ref . . Here we will simultaneously describe the running of the 
scalar masses. For related studies, see also ref. |0. 

2.1 Invitation: the Wess-Zumino Model 

In this subsection we consider a simple example that illustrates many of the main 
ideas we will use in more complicated theories. We consider a massless Wess-Zumino 
model with bare lagrangian 

£0 = jd^ez^^^^+l^jd^e^^^'' + h.c)j, (2.1) 

and higher-derivative regulator terms @| 

Aeg = jd^eZ,^^^^. (2.2) 

We can incorporate soft SUSY breaking by extending the bare couplings A and Zq 
to be ^-dependent (but ^-independent) superfields.[] (A is a superpotential coupling, 
and is not renormalized.) We have regulated the theory in a supersymmetric manner, 
so we can treat the bare couplings as superfields even at the quantum level. 

Because the theory is regulated in a way that preserves SUSY (including the 
spurious SUSY acting on the couplings), the divergences that appear order-by-order 
in perturbation theory can be absorbed by supersymmetric counterterms. That is, 
we can write 

Zo = Z(/i) + 5Z(A,Z(/i),A//i), (2.3) 

where 5Z is the matter wavefunction counterterm. Because the relation between 
the bare and renormalized couplings preserves SUSY, we see that the renormalized 
couplings can also be viewed as SUSY spurious. 

More specifically, we can define the counterterms by computing supergraphs with 
renormalized couplings in the vertices and propagators and choosing the counterterms 
to cancel the divergences. In the SUSY limit where there is no 9 dependence in Zq 
and A, the counterterms have the form [^] 

6C = Jd^eZC{\X\yz^{i2), lAIZ/i)*^*, (2.4) 



^Note that taking a superfield 5' to be x-independent does not violate SUSY, since it amounts to 
imposing the supersymmetric constraint i9^S' = 0. 
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where the form of the function C follows from the fact that theory depends trivially 
on the overall normalization of the fields. 

In the presence of soft SUSY breaking, the renormalized couplings Z and A will 
also depend on 6, and there are new terms in the Feynman rules involving superco- 
variant derivatives acting on the couplings Z and A. However, it is easy to see that 
such terms can be ignored for purposes of computing the counterterms |]rU|. Because 



our regulator preserves the spurion SUSY even in the presence of soft SUSY break- 
ing, we know that the counterterms can still be chosen to be superfield functions of 
A and Z. But local superspace counterterms involving supercovariant derivatives of 
A and Z are forbidden simply by dimensional analysis. We conclude that even in the 
presence of soft SUSY breaking, the counterterms are still given by Eq. ( p.4| ). Note 
what has happened here: the renormalization of the theory with soft SUSY breaking 
is completely determined by a supersymmetric calculation. This is the advantage of 
treating the bare and renormalized couplings as superfields. 

The fact that the theory depends in a trivial way on the scale of the fields can be 
expressed more formally by noting that the the bare lagrangian is invariant under 

$^6^$, Zo e-(^+^')Zo, X^e-'^X, (2.5) 

where A is a 6'-dependent (but x-independent) chiral superfield. The fact that the 
relation between the bare and renormalized parameters preserves this feature can be 
expressed by stating that the renormalized parameter Z transforms the same way as 
Zo: 

Z ^ e-(^+^')Z. (2.6) 

If we view this as a U{1) "gauge" transformation, then \nZ (and In^o) transform as 
a gauge field. This point of view will be extremely useful to us later. 

The relation between the bare and renormalized quantities determined by Eq. (|2T 



Zo = Z(/i) [l + C(|A|Vz3(/i), |A|//.)J , (2.7) 
determines the RG flow of the theory from dZo/dfi = 0. This gives 

= -^^-^C{\X\yz^ |A|/^) = ^{\Xf/Z% (2.8) 

The 6 = 6 = component of 7 is just the supersymmetric anomalous dimension. The 
renormalized soft scalar mass is defined by writing 

Z = z\l- 6WrrA , (2.9) 
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where Z is the renormahzed wavefunction factor. The RG equation for the soft mass 
is determined by the Q'^W' component of Eq. (|2.8| ): 



/^ = - 7(|A|V2^)L.,. = -7'(|A|V^^)^^. (2.10) 
This formula is vahd to all orders in perturbation theory. At the 1-loop level 



1 lA 

^ = "1 

and we recover the familiar result 
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t^f- = T^, AW, (2.12) 

where A = lAjZ"^/^ is the running coupling constant. Eq. (p.8|) also gives the RG 
equation for A terms if we add a non- vanishing 6"^ component to Z{^). 

In the following, we will generalize the procedure followed in this section to general 
renormalizable SUSY theories with soft SUSY breaking. The idea is to include soft 
SUSY breaking by extending the bare couplings Kq to 6'-dependent superfields. As 
long as the theory is regulated in a supersymmetric manner, the bare couplings can be 
viewed as spurion superfields even at the quantum level. We then define renormalized 
couplings K{ij) related to the bare couplings by a superfield relation 

i^o = G(i^(/i),A//i). (2.13) 

The renormalization of the couplings in the SUSY limit then determines the renor- 
malization of the soft SUSY breaking terms as long as the relation does not involve 
supercovariant derivatives acting on K{^). But in a vast class of theories, this is 
guaranteed by simple power counting and symmetry arguments. Eq. ( p. 13 ) therefore 



determines the complete RG flow of all soft SUSY breaking parameters. In the re- 
mainder of this Section, we explain how to carry out these steps for gauge theories, 
which present additional subtleties. 



2.2 Holomorphic Coupling in Supersymmetric QED 

We begin with SUSY QED, a U{1) gauge theory with matter fields $ and $ with 
charges +1 and —1, respectively. This theory can be regulated in a completely su- 
persymmetric manner using a combination of Pauli-Villars fields to regulate matter 
loops and a higher-derivative regulator for the gauge fields. The bare lagrangian can 
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be written as Cq + Aeg, where 

(2.14) 

contains the "physical" couphngs, and 

+ Jd^ew^w^ + h.c., 

G 

contains the regulator terms. Here, Q and Cl are Pauli-Villars fields (odd-statistics 
chiral superfields) and A$ and are cutoffs for the matter and gauge fields, respec- 
tively. We will take the cutoffs to infinity with A$ ~ A^, so there is effectively a single 
cutoff. Note that the bare wavefunction factor Zq appears both in front of the mat- 
ter fields and the Pauli-Villars fields. This is necessary to regularize Zo-dependent 
subdivergences that occur at two loops and beyond. For reference, the components 
of Sq are given by 

where ©o is the (bare) vacuum angle and mxfi is the bare gaugino mass. 

We incorporate explicit soft SUSY breaking by allowing the bare coupling 5*0 and 
Zq to be superfields with nonzero 9 components. Just as in the Wess-Zumino model, 
the fact that the regulator preserves SUSY means that the bare couplings can be 
viewed as superfields at the quantum level, and we can renormalize the theory by 
adding counterterms that are local (in superspace) and gauge invariant. We therefore 
define renormalized superfield couplings S and Z by 

So = S + SS, Zo = Z + 6Z, (2.17) 

where the counterterms SS and SZ are superfield functions of S and Z determined 
order-by-order in perturbation theory to cancel the ultraviolet divergences. 

For Z we can proceed exactly as in the Wess-Zumino model discussed above, but 
we immediately encounter difficulties when we try to renormalize the gauge coupling 
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(2.15) 



(2.16) 



as a superfield. One way to see the problem is that the only manifestly gauge-invariant 
operator that can act as a gauge counterterm is 



C = j(fe\5SW'^W^ + \i.c. (2.18) 



However, the result of a supergraph calculation is necessarily a d^9 integral. At one 
loop, this is not a problem because the one-loop gauge diagrams are independent 
of all couplings (since the gauge coupling is in front of the kinetic term), and the 
counterterm can be proportional to 

Jd^e{D''VWa + h.c.) = Jd^eWWa + h.c. (2.19) 

However, beyond one loop, the coefficient of the counterterm depends on the superfield 
couplings, and the counterterm cannot be written as integral. 

This argument can be sharpened by using the fact that the counterterm 6S is a 
chiral superfield. Because of this, 5S must be a holomorphic function of S, A$, Ac, 
and /i, independent of S"^ as well as Z. We therefore have 

where / is a holomorphic function. Now, the divergence in the gauge coupling g is 
independent of the vacuum angle G to all orders in perturbation theory, since F'^^F^y 
is a total derivative, and therefore irrelevant in perturbation theory.[| Therefore, 

d\m{S) dS ^ ' 

Since / is a holomorphic function, the only possibility is that df/dS is independent 
of S, which implies 

f{S) =a + bS, (2.22) 

where a and b are independent of S. We see that a is the 1-loop contribution, and 
b is identically zero (since the zero coupling limit corresponds to S* ^ +oo). We 
conclude that there is no divergence in the vacuum polarization beyond one loop.0 
If this argument is to be believed, the coupling S satisfies the exact (to all orders in 
perturbation theory) RG equation 

dS 1 , , 

d^ Stt'' 

^We do not address the subtle question of renormalization beyond perturbation theory. 
^Note that this argument does not assume that / is a power series in S. This is important for 
non-abehan gauge theories, where we will see that the perturbation series is non-analytic in S. 
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Figure 1: One-loop diagrams contributing to the vector field propagator. 

This appears paradoxical, since it is known that the /3-function has a (scheme- 
independent) contribution at two loops. 

To understand what is going on, we compute the counterterm explicitly at one 
loop, keeping the couplings as supcrfields. The diagrams are shown in Fig. 1. We 
obtain the contribution to the IPI effective action 

Tipi ^~^J d'^d J d^P V [-f{p^) + SS + SS''] p'^PtV + finite, (2.24) 

where Pt is a transverse superspace projector, and 

d^k -\A^\yz^ -|A$IV2' 



(27r)4 F(A;2 - \A^\yZ^) {k + pfHk + pf - \A^\yZ^) 



1 lA |2/5r2 

In '^^^' + finite. (2.25) 



87r2 -p2 

The IPI effective action can therefore be made finite by adding the counterterm 

55=-^ln— , (2.26) 

where is a renormalization scale. Note that we cannot choose 6S to depend on the 
"kinematic" cutoff |A$|/2, the scale at which the Pauli-Villars regulator cuts off the 
ultraviolet modes, simply because this quantity is not a chiral superfield. On the other 
hand, it is clear that physical quantities depend on A$ only through the combination 
|A$|/Z, together with the bare parameters. This is the key to understanding the 
meaning of the renormalized coupling 5". 

More formally, we note that the bare lagrangian is invariant under 



Zo ^ e-^^+^^^Zo, A$ ^ e-2^A$, So ^ So, 



(2.27) 



where A is a 6'-dependent (but x- independent) chiral superfield. Because Zq measures 
the scale of the fields in the regulated theory, we can choose the subtractions that 
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defined the renormalized Z so that 



Z = ZofiSo + Sl \A^\/Zo, |AG|,/i), (2.28) 

which shows that we can assign the same transformation rule to Z as Zq^ From 
Eq. ( p. 261) we find that the renormalized S{fi) transforms as 



5^ 5- A (2.29) 

Just as in the case of the Wess-Zumino model, we have found a symmetry under 
which Z can be interpreted as a background U{1) gauge field. Eq. ( p.29[) is just a 
reflection of the Konishi anomaly therefore we will refer to this symmetry as 
the (renormalized) "anomalous f/(l)" symmetry. As a consequence of this symmetry, 
physical quantities can depend on 5* only in the combination 



\A^\yz' 



S + S^ -—\nZ = So + Sl + — \n ^ "^2 • (2-30) 



(The right-hand side shows that this combination depends on the kinematic cutoff 
when expressed in terms of the bare parameters.) Notice that 5 — 5*^, which is 
proportional to the vacuum angle, cannot appear in any invariant, consistent with 
the fact that the vacuum angle is not physical in a theory with massless fermions. 
Because of the symmetry defined by Eq. ( |2.27] ) and Eq. ( |2.29D , the relation between 
the bare and renormalized wavefunction factors has the form 

Z„ = ZW/(s + St-^luZ,M^^,lMj, (2.31) 

The RG flow of the theory is determined by dZo/dfi = 0. Due to the loop factor 
multiplying In^ in the above expression, an {n + l)-loop effects are often related to 
n-loop effects. There are many examples of this in the literature, and we also obtain 
new results of this type in subsequent sections. 

Because the correlation functions depend on S only through Eq. ( 2.30 ), the rela- 



tion between the coupling S and a gauge coupling defined directly in terms of IPI 
Green's function is non-analytic in the couplings. As already observed in Refs. , this 
can resolve the apparent contradiction between a holomorphic coupling that runs at 
one loop and the conventional definition of the gauge coupling that runs at all loops. 

We now note that the quantity 

R = S + - ^\nZ (2.32) 



*This may become clearer when we give a IPI definition of Z in the next Subsection. 
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that appears in Eq. ( |2.30|) is a good candidate for a real renormalized superfield 



coupling, i? is a finite quantity tfiat parameterizes the couphngs of the theory, and 
it does not have any unphysical dependence on the scale of the fields. Also, the 
= and 6"^ components of R give the correct RG equations for the gauge coupling 
and gaugino mass to 2- loops. (In fact, Eq. ( |2.32| ) is identical in form to the famous 
equation of Refs. [Q, but note our equation involves only renormalized quantities.) In 
the next Subsection, we will explain the relation between R and a renormalized gauge 
coupling defined from the IPI action, and address the meaning of the 6*^^^ component 
of R. 

We close our discussion of SUSY QED by remarking that there is a completely 
analogous U{1) symmetry with a well defined action on the bare couplings. The 
"gauge transformation" $ ^-^ e"^^ has an anomaly, and so the bare gauge coupling 
must also transform to compensate for the transformation. In our regulator, this can 
be seen from the fact that the Pauli-Villars fields transform under the symmetry, 
so the anomaly can be obtained as the matrix element of the Pauli-Villars mass 
term in a background gauge field. More generally, it is clear that any holomorphic 
regulator yields the anomaly, and the result is that the theory is invariant under the 
transformation 

$^e^$, Zo ^-^ e-(^+^')Zo, So^So-^ (2.33) 
with the regulator Lagrangian invariant. This "bare" or "Wilsonian" anomalous U{1) 



is also a very useful symmetry |12| 



2.3 Real Superfield Coupling in Supersymmetric QED 

We now give another definition of the renormalized gauge coupling, obtained directly 
from the IPI effective action by subtraction at a Euclidean momentum point. This 
corresponds more closely to the "physical" coupling that describes the momentum 
dependence of the effective charge. More to the point, this definition of the gauge 
coupling can be directly understood in terms of component calculations, allowing us 
to make contact between our formalism and conventional calculations. 

In a component calculation, it is natural to define the renormalized gauge coupling 
and gaugino mass in terms of an appropriate IPI correlation function at a Euclidean 
kinematic point. We now show that a definition of this type gives rise to a real 
superfield R whose lowest components are the gauge coupling and gaugino mass. 

Consider the supersymmetric limit first. To define the renormalized gauge cou- 
pling we must consider the gauge invariant bilinears in Wa in the IPI action. Since 
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we include quantum effects we must focus on integrals. By a simple operator 
analysis one finds there exists just one independent term 



TiPi = Jd'p jd^e^{p')w''-^w^ + h.c. + 
= jd^p jd^e li{p^)w"w^ + h.c. + ■ ■ ■ , 



(2.34) 
(2.35) 



where the last identity follows simply by integrating over half of superspace. There- 
fore, 7 can contain the contribution from the tree-level and loop contributions to the 
ordinary gauge kinetic term. We can therefore define the renormalized gauge coupling 
simply by subtracting at a Euclidean momentum point: 



1 



(2.36) 



The role of the operator of Eq. ( p.34|) in generating the all order /3-function was 
already emphasized in Ref. f^. 

We can similarly define a renormalized wavefunction superfield by considering the 
terms in the IPI action that contribute to the matter kinetic term 



ripi= Id^p Id^eCip') $^e^$ + <l>Te'^$ 



and defining 



z = Cip' 



+ 



(2.37) 



(2.38) 



In the presence of soft SUSY-breaking sources in S and Z, the gauge kinetic terms 
in the IPI effective action are 



ipi 



d^p d^9-f{p^)W' 



•ip'^ 



;Wo + \i.c. + 0{Do,S, D^Z,...) 



(2.39) 



where 7(p^) is now a vector superfield function of the couplings 5 + S*"!" and Z, and 
0{DaS, . . .) represents terms involving at least one supercovariant derivative acting 
on the sources. By studying all possible WW and WW terms involving supercovariant 
derivatives, it can be shown that they always lead to terms of second order in the 
soft masses, i.e. they are 0{m? /p^). These terms therefore do not contribute to the 
gauge kinetic term and gaugino mass term in the IPI action. It therefore makes sense 
to define a renormalized superfield coupling by 



i?(/i) = 7(p2 



(2.40) 
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Everything in this definition is manifestly supersymmetric, so the relation between 
this renormalized coupling and the bare couplings is SUSY covariant. The interpre- 
tation of the components of R is given by 

/^'^^^"^^- = [Jd'^O \ (7I0 + 71,2) W'W^ + h.c. 

(2.41) 

The lowest components of R are therefore the coefficients of the gauge kinetic term 
and gaugino mass term, and we identify 

^ =i?(/i)lo, = (2.42) 



Note that this renormalization scheme is mass-independent. 

The O'^O^ component of R multiplies a non-local SUSY-breaking contribution to 
the IPI action. It is instructive to ask what distinguishes this 0{iv?) effect from 
the other 0{m?) WW and WW operators induced by the terms involving covariant 
derivatives acting on the couplings. To do so it is useful to work in components. Since 
there are three component fields A^, A, and D, there are in general three independent 
0{im?) corrections to the corresponding self-energies: 



nT(/) = (pV^-p>^) (1 + ^) 



K 



nA(p')=^(l + ^l (2.43) 




where ka,\,d = 0{m?). A simple operator analysis shows that the terms involving 
super covariant derivatives acting on couplings generate 0{m'^) corrections that always 
satisfy the supertrace sum rule — 4fi;| + = 0. On the other hand the O'^O'^ 
component of R is associated to a non zero supertrace 2 R\q2q2 = 3k^— 4k|+k|,. If one 
computes the effect of the dressed self-energies in Eq. ( p.43|) on the matter self-energy, 
one finds that the only divergent contribution is proportional to the supertrace. This 
simple exercise clarifies why the 9^9^ component of i?, although associated with a 
non-local operator, nonetheless enters into the RG fiow equations of the softly broken 
theory. 
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We now discuss the relation between the real superfield gauge coupling discussed 
here and the holomorphic gauge coupling described in the previous Subsection. Since 
both are perfectly valid parameterizations of the renormalized gauge coupling, we can 
express R in terms of the holomorphic coupling S and Z. The coupling R is clearly 
invariant under the field rescaling Eq. ( p.27| ), so 



= / (^^(/i) + S^f,) - _lnZ(^)j . (2.44) 
Demanding that the holomorphic and real couplings coincide at tree level gives 



i?(/i) = S{f,) + S\f,) - -\ In Z{fi) + ^ + 0{S + S^r\ (2.45) 



where c is a 1-loop scheme- dependent constant. Notice that this expression auto- 
matically gives the correct 2-loop (3 function. Eq. ( 2.45| ) is identical to the famous 
formula of Ref. that relates the IPI and "Wilsonian" gauge couplings. However, 
it is important to remember that the coupling S in our Eq. (p.45|) is a renormalized 
coupling constant. 



2.4 Holomorphic Coupling in Supersymmetric Yang-Mills Theory 

We now consider some additional features that arise in non-abehan gauge theories, 
using the example of a pure SUSY Yang-Mills theory with gauge group SU{N). We 
regulate this theory in a supersymmetric way by embedding it into a finite theory 
with softly broken Af = 2 SUSY. The additional fields in the regulated theory consist 
of a chiral field $ in the adjoint representation (the M = 2 superpartner of the M = 1 
gauge multiplet) and 2N hypermultiplets consisting of chiral fields Vt^ and VLj (J = 
1, . . . , 2N) in the fundamental and antifundamental representations, respectively. 

The bare lagrangian of the theory (written m. M = 1 superspace) is 



(fe Sn tr 



me 



h.c. 



+ ( fd^e V2n'^<^nj + h.c. ) (2.46) 



+ d'e 



+ h.c. 



The coefficient of the Q-^mj interaction is fixed by A/" = 2 SUSY. The TV = 2 SUSY 
is broken explicitly down to TV = 1 by the $ mass term (the mass term for Q"^ and j 
is A/" = 2 invariant). Af = 2 theories are finite beyond one loop ||13[. With our choice 
of matter, the 1-loop beta function vanishes and therefore, in the background gauge. 
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there are no divergences. The parameters and A^ therefore act as cutoffs for 
SUSY Yang-Mills theory, with the fields Q"^ , and Clj playing the role of regulator 
fields. We will eventually take the limit A^, A^ oo with Aq ~ A^, so that there is 
effectively a single cutoff. 

We now show that the finiteness of this theory persists when 5*0 is a chiral super- 
field with nonzero 6 components. Any divergences in the IPI effective action must 
be local (in A/" = 1 superspace) expressions involving the superfield couplings of the 
theory. Because this theory is renormalizable, the divergences must have the same 
form as terms in the lagrangian. There are no divergences when 5*0 is a number, so 
any divergences must be proportional to SUSY-covariant derivatives acting on Sq. 
But such terms have positive mass dimension, so there can be no divergences pro- 
portional to dimension-4 operators. The only remaining possibility is that there are 
divergences proportional to 

Jd^e D^SIQ-^Q J + h.c. or j(feD'^Slii{^'^)+\i.c. (2.47) 

Such divergences can be excluded by considering the (anomaly-free) transformation 

(2.48) 

An ^ A^ ^ e'^^Ac, 

under which D'^Sq is invariant. 

This establishes that the theory above is finite, and therefore provides a regulator 
for the SUSY Yang-Mills theory we want to study. We still need to renormalize the 
theory in order to take the limit Af^, A^ — * oo. The renormalized lagrangian isQ 

£ = Jd^eStT{W"W^)+h.c., (2.49) 

where S is defined by 

So = S + 6S. (2.50) 

The counterterm 6S is fixed order-by-order in perturbation theory to cancel the di- 
vergences as An, Ag oo. 

^The renormalized lagrangian can be thought of as the "effective lagrangian" below the scales 
Aq, Ag- However, we must choose the couplings in the "fundamental lagrangian" Cq as a function of 
Ail and Ag so that the couplings in the "effective lagrangian" are held fixed as the cutoff is removed. 
This can be thought of as "fixing the parameters from low-energy experiment" . 
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At one loop, the vacuum polarization in the background gauge is proportional to 

_ N \MM±^ - in ^ + finite + SS + SSK (2.51) 

where the "physical" cutoff for the $ contribution is |Ag'|/(S' + S"^) due to the non- 
canonical kinetic term for the gauge multiplet. At this order, the theory can be 
renormalized in a holomorphic way by choosing 

^•^ = 7^^^ — + 7^^"^ — ' 2.52 

where /i is a renormalization scale. 

Because the theory is regulated in a supersymmetric manner, the same argument 
used in Sect. 2.2 shows that there are no counterterms beyond one loop to all orders 
in perturbation theory.^ We can therefore choose the counterterm to be given by 
Eq. ( |2.52| ) to all orders in perturbation theory. The renormalized gauge coupling 
defined in this way satisfies the exact RG equation 

dS 3N , , 

As in SUSY QED, the fact that the holomorphic gauge coupling has a 1-loop 
beta function is closely connected to the fact that the subtraction depends on Aq and 
S* + 5*^ separately. Logarithmic divergent loops always involve the "kinematic" cutoff 
|Ag|/(5' + S"^), and therefore the renormalized expansion coefficient is 

s + + ^ .„(S + s.) ^ s„ H- St - m 

-^in^ (2.54) 

We can also define a real superfield coupling from the IPI effective action similarly 
to what was done for SUSY QED. In this scheme, there is a real gauge coupling 
superfield R defined to be the coefficient of the V propagator term in the IPI effective 
action. R must depend on the combination Eq. (|2.54|) , and we find 



R = S + S^ + ^ \n{S + S^) + 0{S + S^)-\ (2.55) 



^Note that the perturbation series is nonanalytic in S, as can be seen from Eq. ( 2.51 ). However, 
the arguments of Sect. 2.2 do not require the perturbation series to be a power series in S, and are 
therefore vahd in this case as well. 
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2.5 General Gauge Theories 



We have so far treated only simple theories where we know how to construct a man- 
ifestly supersymmetric regulator. However, we now argue that our results apply to 
any SUSY gauge theory as long as a supersymmetric regulator exists. The general 
arguments above tell us that the only divergence in the gauge coupling occurs at one 
loop, and has the form 

where Tr is the Dynkin index of the r representation. Here Ac is a cutoff parameter 
for gauge loops and is a cutoff parameter for matter fields in the representation r. 
Note that in order for this formula to make sense, and A,, must be chiral superfield 
spurious, as they are in the examples considered previously. On the other hand, the 
"kinematic" cutoff (the momentum scale at which loop momenta are damped) cannot 
be a chiral superfield, for the simple reason that it must be real. As we have seen, 
Eq. ( p. 561) is consistent with the 2-loop RG equations provided that the kinematic 
cutoff for matter loops is A^^kin = \-Ar\/Zr. The relation between the kinematic gauge 
cutoff and Aq is more complicated, as seen in the example of SUSY Yang-Mills. In 
any case, in order to reproduce the correct 2-loop beta function, physical quantities 
must depend on the combination 

R = S + S^ + ^ \n(S + S^) - Y ^\nZr + 2-loop corrections, (2.57) 

which is the real gauge coupling superfield. In the following we will give further 
evidence for the generality of our conclusions by showing how they arise in dimensional 
reduction, a regulator that can in principle be used for any SUSY theory. 



3 Dimensional Reduction 

So far we have been dealing with regulators that apply only to special theories. How- 
ever, in order to be able to calculate higher order effects in any theory, including 
the supersymmetric extension of the Standard Model, the only practical regulator is 



dimensional reduction (DRED) |jT^ . In this section we show how the holomorphic 
and real gauge couplings arise in DRED. We also show that the procedure of analyt- 
ically continuing the renormalized couplings into superspace picks out the so-called 



DR scheme [0 in which the e-scalar mass does not appear in physical quantities. 
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3.1 Real and Holomorphic Gauge Coupling in Dimensional Reduction 

The renormalization of SUSY gauge theories in the framework of DRED was clarified 
more than a decade ago by Grisaru, Milewski and Zanon (GMZ) They pointed 
out that in d = 4 — 2e dimensions, there is an additional supersymmetric and gauge- 
invariant local operator 

OGMz = ^?rtr(r^r,), (3.1) 

where g^"^ is the metric in the 2e "compactified" dimensions, and is the superfield 
gauge connection defined by 



= (e-^Z^"e^) . (3.2) 

This operator is an 0{e) (or "evanescent") operator, with the property that 

J d^e Ogmz = ejd^e tr (1^"1^„) + h.c. (3.3) 

(Note that g^Tfj^Ti, is real.) Therefore, the quantity / d'^O Cqmz is a dimension-4 term 
that can appear as a counterterm for the gauge kinetic term. 

Taking this into account, the bare lagrangian is 

Co = Qd'^e SotiiWWa) +h.c}j + y"d^^To^f^tr(r^r^) + matter terms. (3.4) 

We can incorporate soft SUSY breaking by extending 5*0 and Tq to 6'-dependent 
superfield spurious. Because DRED preserves SUSY, we can treat So and Tq as 
superfields even at the quantum level. The meaning of the higher components of Tq 
is given by 



Jd'eTogrtT{T,r,) = e \^Jd^e (ToI + To\,,) tiiW^^W^) + h.c. 



(3.5) 



That is, the lowest components of Tq are contributions to the gauge coupling and 
gaugino mass, and the O'^O'^ component is the e-scalar mass. 

We now renormalize the theory by writing 

So = ^i-^'{S + 5S), To = /i-2^(r + (5T), (3.6) 

where 5S and 5T are counterterms that are determined order by order in perturbation 
theory to absorb the 1/e divergences. Note that we include a finite renormalized 
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value for T. This corresponds to including evanescent effects: the scalar and Q'^ 
components of T are C(e) contributions to the gauge coupling and gaugino mass, and 
the (P'Q'^ component of T is a renormalized e-scalar mass parameter. We will return 
to the significance of these parameters below. If we compute using supergraphs, all 
divergences appear in the IPI effective action in the form / d'^OO/e^, where (9 is a 
local (in superspace) gauge-invariant supersymmetric operator, so the counterterms 
can be defined to preserve the SUSY acting on the coupling constants. 



Ref. 10 show that at one loop, the divergences can be absorbed in 5S, but at two 
loops and higher, all divergences must be absorbed in 5T. This sheds considerable 
light on the origin of the 2-loop running of the gauge coupling, as follows. At 2 loops 
(and higher), a pole in 5T will appear as a result of subdivergences. By Eq. ( |3.3|) , 
this corresponds to a 1/e pole in the counterterm for the gauge coupling, which affects 
the beta function. The fact that a pole arises only from subdivergences explains 
why the higher-loop contributions to the gauge coupling beta function are determined 
by the anomalous dimensions of the matter fields. 

New features arise if we include soft SUSY breaking by extending the couplings to 
superfields. At one loop, we find an ultraviolet divergent contribution to the e-scalar 
mass: 



(3.7) 



Although this is a finite effect, it is known that renormalization of the e-scalar inter- 
actions is required to preserve unitarity |jl8|, |l9l. (Indeed an explicit calculation of 
Poppitz and Trivedi |^ shows that infrared divergences arise at 2-loops if the e-scalar 
mass is not renormalized.) 

To subtract the divergence in the e-scalar mass in a way that preserves SUSY 
acting on the couplings, we must add the 1-loop counterterm 



TT^ e 



(3.8) 



The logs ensure that the counterterm for the e-scalar mass has the correct dependence 
on the gauge coupling and is independent of the wavefunction of the matter fields. 
Note that the scalar and 6^ components of 6T give rise to finite contributions to the 
gauge coupling and gaugino mass. This restores the dependence of the renormalized 
gauge coupling on \nZ and ^(5* + 5^). 

We now have all the ingredients we need to define the renormalized holomorphic 
and real gauge coupling superfields in DRED. The holomorphic gauge coupling is 
defined simply by S. Because 6S contains only 1-loop divergences (and 5*0 is n 
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independent), S runs only at one loop. On the other hand, because of the subtraction 
in Eq. (|3.8| ), the components of S do not give the renormalized gauge coupling and 
gaugino mass. Rather, these are given by the lowest components of a superfield R, 
defined by 

R= S + + eT + 6T^^\ (3.9) 

where ST^^^ is the coefficient of 1/e in 5T. From Eq. ( |3.8| ), we see that the quantities R 
and S satisfy precisely the relation derived in the previous Section for other regulators 
and renormalization schemes: 

R = S + S^ + ^ \n{S + S^)-J2^\nZr + + S^y'). (3.10) 



The definition Eq. (|3.9|) also shows that physical quantities must depend on S through 
R, since it is the components of R that multiply the kinetic terms and gaugino mass 
terms in the lagrangian. 

We need to understand what scheme in component calculations is picked out by 
the procedure above. It is useful to define a bare gauge coupling superfield 

R, = So + Sl + eTo (3.11) 

in terms of which the bare gauge coupling and gaugino mass are 

~ = -^olo 5 J~ = i^O + ^TQ)\g2 = Ro\g2 ■ (3.12) 

while the renormalized couplings are (see Eq. ( p.9| )) 



— -R\q, -—^-R\g2. (3.13) 

The relation between the bare and renormalized couplings is therefore determined by 
the components of 

fio = f.-^'(ii+^ + i:^). (3.14) 

where JT*^") is the coefficient of l/e"" in ST. We assume that SS and 5T consist 
of pure l/e poles. This corresponds to modified minimal subtraction (MS) if we 
rescale /x appropriately, writing fi = fl^Je'^/An and writing all expressions in terms of 
jl. Eqs. ( |3.12| ) and ( p.l3D then show that g and mx are precisely the renormalized 
couplings in DR. 
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When we consider the inclusion of matter with soft scalar masses, the scheme 



picked out by the procedure above is identical to DR |jT6[. To understand the issues 
involved, note that there appears to be an extra renormalized parameter in DRED, 
corresponding to an e-scalar mass. This parameter has a non-trivial RG evolution, 
and so cannot be set to zero at all scales. However, the e-scalar mass is an evanescent 
effect, and does not give rise to an additional parameter at the quantum level. The 
way this works is that if we renormalize the theory with an arbitrary e-scalar mass 
parameter, it only appears in physical quantities in the combination |l^ 



(3.15) 



One can then define the scheme DR by declaring the combination above to be the 
renormahzed soft scalar mass. DR' is therefore the scheme in which the e-scalar mass 
does not appear in any renormalized expression for arbitrary values of fi. 

In terms of the superfield couplings, the renormalized e-scalar mass corresponds 
to the term eO'^O'^ in R. But because we subtract all the 1/e poles in R, the IPI action 
is a finite function of R. Therefore, there is no explicit dependence on m\ in physical 
quantities, for any value of fi. This is sufficient to prove that the scheme we have 
defined is identical with DR'. Our procedure extends the definition of DR', given in 
ref. [|16] at the 2-loop level, to all orders in perturbation theory. 

Note that the inclusion of the evanescent eT term in (|3.9| ) is essential for R to 
satisfy the d-dimensional RG equation 



fi^ = 2eR + (3{R). 



(3.16) 



This is easy to check at lowest order by considering the RG equation for T. Therefore, 
in our scheme plays a role similar to that of the 0{e) term in the ci-dimensional 
RG equation for g'^: it insures fi independence of the bare coupling Qq, but is irrelevant 
in calculations. 



To see more explicitly the connection to naive DR, consider the relation between 
the bare and renormalized wavefunctions for the matter fields 



■'r,0 



1 + E 



5Z^"\R) 



n=l 



(3.17) 



Taking the 9 9 components of both sides gives 



r,0 



— 



d 

dR 



5Zl^\R) 



rrij^ + 1/e poles. 



(3.1^ 
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In our scheme, the renormahzed scalar mass is ml = — lnZ\g2(j2, while the finite term 
on the right-hand side is the scalar mass in DR {not DR'), since it corresponds to 
minimal subtraction. Comparing Eqs. (|3.18|) and ( |3.15|) , we see that mf. is identical to 
m^— / to 2 loops. (But note that our scheme is defined to all orders in perturbation 
theory.) 

Let us summarize the main results. In the supersymmetric limit where the explicit 
soft breaking is turned off, we can renormalize the theory by (modified) minimal 
subtraction, defining renormalized couplings in the DR scheme. Our result is that if 
we include renormalized soft terms by analytically continuing both the renormalized 
couplings and the counterterms (defined as functions of the renormalized couplings) 
into superspace via 

\^R, Zr ^ Zr, (3.19) 

this defines a valid subtraction scheme for the softly-broken theory. This picks out a 
unique scheme for the soft terms to all orders in perturbation theory, which we call 
SDR for supersymmetric dimensional reduction. (At two loops, SDR coincides with 
DR', so we can think of it as an all-orders definition of DR'.) In SDR, the RG equa- 
tions for all soft parameters is determined by the RG equations in the SUSY limit, 
to all orders in perturbation theory. For instance, in gauge mediated models (see the 
next section), the analytic continuation of Eq. ( |3.19|) is simply performed by substi- 
tuting M — > M + in the effective couplings of the low-energy supersymmetric 
Standard Model. 

We close with two comments on the superfield coupling R defined above. Note 
that the finite 6^6^ component of R defined in DRED corresponds to an infinite 
contribution to the e-scalar mass. In our definition of R from the IPI effective action, 
the 9"^ 9"^ component of R was related to a nonlocal effect. It is interesting to see the 
connection between these effects explicitly by considering softly broken SQED as in 
Section 2.2, but dimensionally reduced to 4 — 2e dimensions. After subtracting the 
Z independent 1/e divergence the gauge self-energy has the form 



ripi = — ld^9 \nZ 



-g^" tr(r^r, + h.c.) + tr W^W^ + h.c. 



(3.20) 



-|- (Z-independent) + 0{DaZ, 



If we write this out in terms of components of Z, we see that the the terms involving 
Z\ and Z\g2 are local and exactly cancel between the two terms in brackets. What 
is left, from \nZ\g202, is just a divergent e-scalar mass, see Eq. ( |3.7| ), and a non- 
local correction to the gaugino self-energy, see (|2.43|) . Anyway, we must subtract the 
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divergent e-scalar with a superfield counterterm as (|3.8| ), so that in the subtracted 
IPI, the dependence on In Z is all coming from the non-local operator. This shows 
that the "chiral" components of R defined in DRED and by IPI subtraction differ 
only by finite analytic (^^-independent) terms, that is, by a change in scheme. In this 
sense, the two definitions are equivalent. 

A closely-related issue involves the relation between the origin of the XwZ term 
in R in DRED and in the general discussion of SQED given earlier, where it was 
inferred from the anomalous f/(l) symmetry. It is conventionally said that there is 
no rescaling or chiral anomaly in DRED, and it may appear that there is no direct 
connection between these arguments. However, an intriguing clue can be seen by 
considering the bare Lagrangian with couplings 5*0, To, and Z^. This Lagrangian has 
the symmetry 

To^To + A + A+, 5oH^5o + eA, Zo^Zo, (3.21) 

which ensures that physical quantities depend on the combination So + 5*0 + eTo. 
However, arbitrary values of To lead to inconsistencies (loss of unitarity and IR di- 
vergences). Up to two loops the choice 

To = --\-lnZo (3.22) 

eliminates the problems. But with this choice, physical quantities depend on the 
combination 5*0 + — In 2o/47r^, which is just what is required to obtain the anoma- 
lous f/(l). We believe that these are very suggestive connections that come close to 
exposing the anomaly in DRED, and we plan on exploring this point more completely 
elsewhere. 



3.2 Two-loop Renormalization Group equations in DR 



We can check explicitly that the scheme defined above is equivalent to DR at NLO 
by computing the 2-loop RG equations for the gluino and sfermion masses. Consider 
the real gauge coupling, given by 



R{^^) = sif,) + ^t(^) + 1^ 



In 



^(/i) + 5t(/.) 



T 



InZrifi), 



(3.23) 



where S is the holomorphic gauge coupling. The gaugino mass is given by m\ 
\nR\g2, so its NLO /? function is easily derived from Eq. ( p.23| ): 



/i- 



in 



2\2 



TGb-2Y,TrCr]mx. 



(3.24) 
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where b = 3Tg — I]r ^r- This equation agrees with the exphcit component calculations 
in DR. A similar derivation, based on the Konishi anomaly, was given by Hisano and 
Shifman 0. A new feature of the present treatment is that R also governs the 
evolution of the dimension-2 soft terms. To see this, consider 

1 



R 



(3.25) 



According to our discussion above, -R|e202 corresponds to a 1/e counterterm for the 
e-scalar mass. Eq. ( |3.25 ) agrees with what is found in explicit component calculations 
(Notice that the quantity on the right-hand side is proportional to the supertrace 



weighted by the Dynkin indices.) Now, consider the 2- loop RG equation for matter 



fields in DR m 



dlnZr 



/i- 



djj, 



2Crg^ + [3Tg - T 




(3.26) 



where T = J2r Tr- Its continuation into superspace simply amounts to the substitution 
g"^ ^ 1/R, Z —>■ Z. The RG equation for the scalar masses is then obtained by taking 
the 6^6^ component of Eq. (|3.26|) . This gives 



drn^ 



87r2 



+ 6(3Tg 



(3.27) 



2Cr)ml 



which agrees with the result in DR [||, |l6l. The same check can be done for 

the evolution of A- and i?-terms and in the presence of Yukawa interactions. 



4 Gauge-mediated Super symmetry Breaking 



We now show how to apply the formalism of the previous Section to perform cal- 
culations in gauge-mediated SUSY breaking (GMSB) models. We begin by briefly 
reviewing the calculation of the leading gaugino and scalar masses in GMSB, as per- 
formed in Ref . |^ . We then turn to new calculations at higher loop orders. The main 
new result in this Section is that the gaugino masses are insensitive to the couplings 
in the messenger sector up to four loops. This "screening theorem" means that it is 
possible to make precise predictions for gaugino masses even when the SUSY break- 
ing dynamics is strongly coupled. The scalar masses are not screened in this way, 
and are therefore sensitive to strong SUSY-breaking dynamics. We also compute the 
NLO corrections to SUSY-breaking masses in GMSB, which correspond to 2-loop 
corrections for gaugino masses and 3-loop corrections to the scalar masses. 
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4-1 Leading Results 

In this Subsection, we briefly review the main results of Ref. for completeness. 
Consider the fundamental theory 



c = d^e 



+ Jd^eS'tT{W''Wa)+h.c. (4-1) 
+ Jd^e XXQQ + h.c., 



where Q,Q are the messengers, are observable sector fields, and X is a singlet. X 
is a background chiral superfield that parameterizes the effect of SUSY breaking via 

XX = M + e^F, (4.2) 

with the assumption F -C M^. Our notation is appropriate to the case where there is 
a single gauge group, but our formulas are trivial to generalize to the case of product 
gauge groups. Below the scale M, the effective lagrangian is 

c = jd'eY.Zrqy''\, 

(4.3) 

+ jd^eS ti^WW^) + h.c. + • • • , 

where the omitted terms consist of higher-dimension operators. The low-energy gauge 
coupling is given by tree-level matching and one-loop running to be 

S{^^) = S'if^o) + t|-^ 1^ - + TTT^ 1^ (4-4) 

where 

b' = b-N, b = 3TG-Y.Tr, (4.5) 

r 

are the beta function coefficients in the full and effective theories, respectively. = 
J2q Tq is the "messenger index" . Here fiQ is an ultraviolet scale where the theory is 
defined; this means that we must evaluate derivatives holding the running couplings 
at the scale /iq fixed. 

The dependence of the low-energy effective Lagrangian on the SUSY-breaking 
effects is given simply by making the replacement 

M ^X (4.6) 
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in the dependence of the effective couphngs S and Zr. (Notice that to simphfy the 
notation we have absorbed A in the definition of X) . It is this "analytic continuation" 
that is at the heart of the method of Ref . ^ . We can now read off the gaugino mass 
from 



dX 



16n^ M' 



(4.7) 



where the notation "|o" denotes setting 9 = 9 = and X = M. Note that this 
automatically gives the correct RG improvement of the gaugino mass. Eq. (|4.7] ) 
involves the holomorphic gauge coupling, which is equivalent to the real superfield 
coupling at one loop. The use of the real gauge coupling is crucial for the higher-order 
calculations we do later. 

We now consider the contribution to the gaugino mass coming from higher- 
dimension operators in the effective lagrangian 0. Operators in the effective la- 
grangian consist of analytic terms in the light fields and the background X and their 
derivatives divided by powers of X. The lowest-dimension operator respecting the 
U{l)j^ symmetry that can contribute to the gaugino mass is 

'XW^X 



6C 



eg 
167r2 



d''9 



\X\ 



tr(iy"W„) + h.c. 



Eq. 



]8D gives a contribution to the gaugino mass of order 



5m \ 



mx 



\M\' 



(4. 



(4.9) 



This is negligible if F <^ M^. It is easy to see that all other higher-dimension 
operators also give contributions to the gaugino and scalar masses that are suppressed 
by powers of F|-/M^. 

We now turn to the calculation of the scalar mass, where the correct continuation 
into superspace is M — >^ VXJO 0. We compute the matter field wavefunction 
coefficients Zr, whose ^-dependence contains SUSY breaking from the dependence on 
the threshold at M: 



mlifi) = - 



(9^ lnZr{fi) 



dXWX 



1 dHnZrifi) 
1 {dln\X\y 



M 



(4.10) 



where we have used the fact that In^,, is a vector superfield, and therefore depends 
on X through The 1-loop RG equation for Z^ is 



7r = Z^- 



dlnZr{fJ.) Cr 



47r2 S + S^ 



Ir 



rflnZ;(/i) Cr 



47r2 S' + ^'t ■ 



(4.11) 
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Computing using 1-loop running and tree-level matching, we have 



(4.12) 



This gives 



91nZr(/i) /■^^ d[i! 97r.(/i' 



ainixi 



|x| /i' 91n|X| 
Cr dfi' d 



47r2y|x| /i' 91n|X| V^(/i') + -^K/^O 



(4.13) 



Note that the explicit |X| dependence from the limits of integration cancels in the 
derivative because of the tree-level matching conditions. From Eq. we see that 



(4.14) 



h' X^X h /y^ 



167r2 /ig 167r2 XtX' 
which depends on X only through |X|, as required. We then obtain 

dlnZrijj) Cr ff' dfi' b' - b f 1 



dln\X\ An^Jx fi' Svr^ \S{fi') + S^fi') 

Computing one more derivative yields 



(4.15) 



9^ lnZr(/i) 



a(ln|X|)2 



2ax 



(Stt 



2^2 



(4.16) 



At=M 

where we used the definition of the messenger index N = b — b'. This gives a scalar 
mass 



2,.., CrNa^{M) 



mt{M) 



M 



(4.17) 



It is remarkable that the finite part of a 2-loop graph can be evaluated from a 1-loop 
RG computation. In the present approach, this arises because Zr depends on |X| 
only through the values of running couplings, and derivatives with respect to |X| 
therefore bring in extra loop factors. 



4-2 Gaugino Screening 

We now consider corrections to the gaugino mass. Very generally, we will find that 
contributions from messenger interactions to the gaugino mass are suppressed by 
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additional loop factors beyond the naive expectation, a result we refer to as "gaugino 
screening". We will see in Sect. 4.4 that the scalar masses are not similarly screened. 

The main point is that the holomorphic gauge coupling is given exactly by 

b' -b 

S(i_i) = ^InX + (X-independent), (4.18) 

where b {b') is the beta function coefficient in the effective theory below (above) the 
messenger scale. (If the SM gauge group has a standard embedding into a larger 
messenger group above the messenger scale, then b' is the beta function of the larger 
group.) The physical gaugino mass must be read off from the ^-dependent compo- 
nents of the real superfield gauge couphng. (As explained above, the holomorphic 
gauge coupling has unphysical field rescaling invariance that is not present in physi- 
cal quantities.) The real gauge coupling is related to the holomorphic gauge coupling 
by 



i?(/i) = ^(^) + 5t(/i) + ^ln S{f,) + S\f,)\-Y.^lnZ, 



T 

J. r 



87r2 



(4.19) 



+ OiS + S^)~\ 



The dependence on the wavefunction factors Zr contains the information about the 2- 
loop RG behavior of the physical couplings. Since S is just given exactly by Eq. ( |4.18| ), 
and since the sum over r runs only over the light fields, R is not affected at the NLO 
by the messenger interactions. That's all there is to the proof! 

Because the leading dependence on the messenger interactions comes from Zr in 
Eq. ( [4.19[ ), it is easy to see that 



5mx f g 



r\-/ 



mx \47r 



4 



g2 M'' 

^mess I 1 

16^ IT 



(4.20) 



The {g/Air)^ factor arises because the messenger fields interact with matter only at 
2 loops. The first term in square brackets represents a threshold correction due to 
a messenger coupling (yfmoss, "while the term ln(M'/M) represents the sensitivity to 
mass splittings among the messengers. Such mass splittings will arise if the various 
messengers have different Yukawa couplings A to the same source X (see Eq. ( [1.1[ )). 
In the next Subsection, we perform explicit calculations of the gaugino masses and 
NLO, and we will see how the screening theorem manifests itself in detail. In the 
remainder of this Subsection, we confine ourselves to some qualitative remarks. 

Consider for example the dependence on the messenger Yukawa coupling A. At 
leading order, the low-energy gaugino masses are independent of A, but one may 
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naively expect important quantum corrections if A is large. This is not an artificial 
possibility: if the Yukawa coupling arises from composite dynamics, the value of A 
will be close to the perturbative limit A ~ 47r at the compositeness scale |2^. In 
this case, g^^^^/ {16n^) ~ 1 in Eq. ( [4.20| ), but 6mx/mx is still suppressed by two weak 
loops. Therefore, the gauge-mediation gaugino mass relation are rather insensitive to 
strong dynamics of the messenger fields even if A is close to the perturbative limit. 

Another interesting example is the case in which different messenger fields have 
different Yukawa couplings to the same supersymmetry breaking source X. In other 
words, the various messengers have different masses M but the same ratio F/M. For 
example, in a GUT model with a messenger scale much lower than the GUT scale, the 
running of the messenger Yukawa couplings between the GUT scale and the messenger 
scale can induce splittings of the messenger masses of order ((/^/(IGtt^) In^McuT/M), 
which can be 0{1) even if the messenger Yukawa interactions are unified at the 
GUT scale. Now, Eq. ( [4. 201 ) shows that, even for 0{1) messenger mass splittings, 
the minimal GMSB relation between the different gaugino masses are only violated 
by 0{{g/4:7i)^). Therefore, the gaugino masses do not depend on the assumption of 
universality of the messenger Yukawa couplings at the messenger scale even at NLO, 
as long as the Yukawa couplings are of the same order and there is a single source X 
of SUSY breaking. 

Similar considerations apply to models with vector messengers. In such models, 
the vacuum expectation value that breaks SUSY also breaks a larger gauge group 
down to the standard-model subgroup. There are therefore massive gauge bosons 
charged under the standard-model gauge group that act as SUSY-breaking messen- 
gers. Ref. computed the leading contribution of vector messengers to the scalar 
and gaugino masses, and showed that the contribution to the scalar mass-squared is 
negative. The leading contribution to the gaugino mass from the vector messengers 
also arises at 4 loops, and again has the order of magnitude given in Eq. (^4.20|) , where 
(7mess is now the messenger gauge coupling. This is important because the messenger 
gauge coupling can be strong at the messenger scale. (For example, this occurs in the 
models of Ref s. pi). 



4-3 Gaugino Masses at the Next-to-Leading Order 

We now compute the NLO corrections to the gaugino masses in DR. In compo- 
nents these corrections correspond to threshold effects at the messenger mass scale 
described by 2-loop Feynman graphs, together with the two-loop RG evolution from 
the messenger scale to the physical scale. In our approach, these corrections can 
be extracted from the expression of the real superfield R. In DR, the NLO match- 
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ing at the messenger scale is simply obtained by requiring continuity of -R(yu) at the 
threshold of the physical messenger mass ||2^ 

Here Zm is the wavefunction factor for the messenger fields. Following the notation 
introduced in sect. 4.1, primed (unprimed) quantities refer to the theory above (below) 
the messenger mass scale. In terms of the value of R' at an arbitrary high-energy 
scale fj,o, much larger than the messenger scale fix, at the low-energy scale fi we find 

R(^) = R M + 16^ 1- ;^ + 8^ 1- ^^siiI7) 'V^' zAj^y ^ ^ ^ 



167r2 /ig 87r2 Re^'(;Uo) 

Here is the gauge coupling at one loop (see Eq. ( [4.14| )), and R'{fio) = R-e5"(/io) 
gives a SUSY-preserving boundary condition on the gauge coupling. The sums in 
the previous equations extend over the different matter superfields. Substituting 
Eq. (|;22D into Eq. ( ^:23D , we obtain 

i?(/i) = i?'(/io) + -TT^ In ^ + -TT^ In ■ 



167r2 /i2 I67r2 filZlj{fio) 

^ 8^ Re^'(/io) ~ ? 8^ Z;(/io) ■ ^^'^^^ 

Notice that in this expression the explicit dependence on Zm{iJ'x) has dropped out. 
An implicit dependence appears from higher-order contributions in the matter wave- 
function renormalization Zr{fi). However, the NLO expression for the gaugino mass, 
which requires only the leading contribution to Zr{iJ,), is independent of ZM^fix)- 
This is a manifestation of the "gaugino screening" theorem discussed in Sect. [4.2| . 
See see that at this order in perturbation theory, the gaugino masses are not affected 
by new messenger interactions. Similarly, the VT-ino and 5-ino masses have no 0:3 
corrections from messenger thresholds, but only from their RG evolution below the 
messenger mass. 

To obtain the expression of the gaugino mass, we take the F component of 
Eq. i ^Ml - 
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l-g\^^)TG/{S^l') 

This equation gives the NLO expression of the gaugino mass in terms of the SUSY- 
breaking part of the hght matter wave functions Zr{fi) at the leading order. To 
complete the calculation, we now compute In .2,. (/i) 1^2 for matter fields including both 
gauge and Yukawa interactions. For simplicity, we give the result for a simple gauge 
group, but the generalization to a product group is completely straightforward. The 
relevant 1-loop RG equations are: 
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dr 2 



(4.26) 
(4.27) 
(4.28) 



where y is the running Yukawa coupling (physically normalized by appropriate wave- 
function factors). Here, dr is the number of fields circulating in the Yukawa loop, 
and 



(4.29) 



with the sum extended to the field participating in the Yukawa interaction. If g is 
the QCD coupling, and y is the top-quark Yukawa coupling, we have C = 8/3 and 
D = Q. Taking the F component of the solution of Eq. ( |4.26| ) for Zrln), we obtain 
the final expression for the gaugino mass including QCD (as) and top-quark Yukawa 
{at = yl/{'iT[)) corrections 
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I{i)Y.Trdr 



(4.30) 



(4.31) 



where the sum is taken over the colored light fields. This result has been recently 



confirmed by an explicit component calculation [^. Notice that in the above equa- 
tions the dependence on the physical messenger mass appears via ^, and it is of the 



form given in Eq. ( 4.20 ). 
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In order to obtain the pole gaugino mass we have to include also the finite one- 
loop corrections at the infrared threshold. For the gluino, in the DR scheme they are 
given by ^7 
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(4.32) 



A3. 



J^(x) = l + 2x + 2x(2-x)lnx + 2(l-a;)Mn|l-a;| 



(4.33) 



The function includes the effect of the gluon-gluino and quark-squark loops in the 
approximation in which all squarks have equal mass rhq. Since we have neglected 
weak corrections, the SU{2) x f/(l) gaugino masses receive no contributions from 
infrared thresholds. The final expressions for the three gaugino masses improved by 
as and at corrections are then given by 
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(4.34) 
(4.35) 
(4.36) 



where 02 = foi at the unification scale. 
The NLO correction to the gluino mass 

\LO 



pole \ NLO 

A3 ) 



/("^Ag 



\LO 



1 is shown in Fig. 



We have assumed {m\^ = 600 GeV and tan P = 2, but the result is very insensitive 
to this choice. In particular, the value of tan P is unimportant because the top- 
quark Yukawa contribution in Eq. ( [4.34| ) is negligible. The NLO contribution from 
messenger loops, which is obtained by setting ^ = 1 in Eq. ( |4.34| ), is about +A- 
5%. However, the NLO gauge RG evolution contributes a negative contribution (see 
Eq. ( [4.34[ ) and Fig. 2) that almost completely cancels the messenger contribution for 
very large running (M ~ 10^^ GeV). The finite gluon-gluino loop gives also a large 
positive contribution of about +10-12% to the gluino mass. This effect is partially 
compensated by the quark-squark loops, if the ratio m?/M| is not large, as in the 
case of several messenger flavors {N > 1). This explains why the NLO correction to 
the gluino mass is very important for small M and N, but signiflcantly decreases for 
larger values of M and N (see Fig. 2). 
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Figure 2: NLO correction to the gluino pole mass, as a function of the messenger 
mass scale M , for messenger index = 1, 3, 5. We have taken a leading-order value 
of the gluino mass of 600 GeV and tan/? = 2, hut the results are rather insensitive 
to these choices. The curves are interrupted at values of M that require F = to 
obtain the required gluino mass. 
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The QCD corrections to the SU (2) x f/(l) gaugino masses vanish at the messenger 
scale, as expected from the "screening theorem" previously discussed. The effects from 
the RG running, shown in Eqs. ( [4.35| )-( p:.36| ), tend to cancel between the gauge and 
Yukawa term, and give a contribution to the weak gaugino masses that is at most of 
few percent. 



4-4 Scalar Masses at the Next-to- Leading Order 



We can now also compute the NLO corrections to the squark and slepton masses in 
DR', which correspond to 3-loop diagrams. The RG equation for the wave-function 
renormalization of a matter field r is 

d 



dfi 



In Zr 



(4.37) 



The gauge contribution to the anomalous dimension 7^ at the NLO is given by ^T|, ^ 



7r 
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Cr- — 5- + Cr [3Tg — 2Cr — T] 
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The SUSY-breaking scalar mass is obtained from Eq. ( [4. 101) 

IdHnZrin) F ^ 



(4.38) 
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{d\n\X\y 



(4.39) 



where 'jr il'r) is the anomalous dimensions below (above) the physical messenger scale 
fix (see Eq. ( [4.21| )). Note that 7.^ in the low-energy theory depends implicitly on fix 
from the matching conditions at the messenger threshold. Notice also that the lowest 
matching correction for the wave function at the messenger scale fix is at 2-loops. 
This corresponds to the addition of an 0{a{X)'^ /IGir"^) term inside square brackets 
in Eq. (|4.39| ). The resulting correction to the squark mass is (9(a^). 

For simplicity, we will give the expression of the scalar masses evaluated at the 
messenger scale, as the 2-loop running from fix to the low-energy scale fi is well 
known p2|, [10|, |l^]. In this case, the action of d'^/{dln |X|)^ on Eq. ([4.39|) gives, at 



the NLO in gauge interactions, 
F 2 dlnfix 



m^{fi 



M 
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d\n\X\ dlnfix 



H{B'{fix))-lr{R{fix))-lr{R{fi))] 
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(4.40) 



Here 7m = In /c? In /i is the anomalous dimension of the messenger superfield 
at the leading order, which depends not only on gauge interactions, but also on any 
new additional interactions of the messengers. In particular, including the Yukawa 
interaction in Eq. ( |4.1| ), we find 
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This explicitly shows that the "screening theorem", valid for gaugino masses, does 
not apply to scalar masses. 

We can now evaluate the derivatives of R, using the expressions obtained in the 
previous section: 
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(4.43) 



Notice that in Eq. ( ^.42| ) we have kept only the leading term in the perturbative 
expansion, since in Eq. (|4.4CI| ) it multiplies the factor 9(7^ — 'jr)/dg^, which is an 
NLO quantity. Putting it all together, we obtain the final expression for the scalar 
masses at the NLO 
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(4.44) 



Assuming that the messengers belong to fundamentals of SU{5), the NLO expres- 
sion for the QCD contribution to squark masses is 
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(4.45) 



Here A3 is the messenger Yukawa coupling for the color triplet. Notice that 0{al) 
contribution to squark masses from the messengers tends to cancel the contribution 
from gauge and matter fields, as long as N is not too large. NLO corrections to 
slepton masses from QCD and new messenger interactions come only from the factor 
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(1 — 7m) in Eq. ( [4 .401) . Since, in our case, weak-doublet messengers are color neutral, 
the SU{2) contribution to left-handed slepton masses is corrected only by the factor 
(1 + A2(/ix)/87r^). However, for a generic choice of messengers, the QCD corrections 
are non- vanishing. Notice also that in general A2 7^ A3, although they may be related 
in a GUT. Finally, the improved expression for the right-handed slepton mass is 
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(4.46) 



In Eq. ( [4.46| ), fix can correspond to the mass scale of either the triplet or the dou- 
blet messenger mass. The difference between the two definitions is 0{al), which is 
negligible in our approximation.!] On the other hand, fix in Eq- ( [4.45|) has to be 
interpreted as the triplet messenger mass, since we include terms (9(a|). 

In conclusion, because of the absence of a "screening theorem", the NLO cor- 
rections to scalar masses are quite dependent on the model assumptions. They are 
sensitive to new messenger interactions, like the messenger Yukawa couplings, and 
they depend on the messenger representations in a way that cannot be described only 
by the messenger index N. 



4-5 D-type Supersymmetry Breaking 

We now consider leading SUSY breaking effects in theories where the dominant source 
of SUSY breaking is a D-type soft mass for the messengers rather than a F-type 
mass, as considered previously. Some of these results have already been derived 
in the language of renormalized couplings in Sect. 3.2. We discuss them here in a 
manifestly "Wilsonian" picture, that is, by simply computing in the theory with given 
bare parameters. We do this in part for variety, and in part to show how these results 
follow from the "Wilsonian" anomalous U{1) symmetry.^ Consider a gauge theory 
with bare lagrangian 

£0 = jd^eSQii{W''W^)+h.c. + yc/^^Z,,o$Je^"''$r, (4.47) 

regulated in a supersymmetric manner. Assume that the theory contains bare soft 
masses, parameterized by 

-eWml^. (4.48) 



2, 



r,0 — ^r,0 



^Higher orders in the electroweak couplings can be computed following the same procedure used 



to obtain Eq. (4.44), with the introduction of separate messenger thresholds. For an application of 
the method of Ref. |^ to the case of multiple messenger thresholds, see Ref. [ p8[ . 

*This symmetry is extremely useful in obtaining physically interesting results for non-holomorphic 
soft terms in strongly coupled SUSY gauge theories with small soft breakings S . 
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As discussed above, this theory is invariant under the "Wilsonian" anomalous U{1) 
transformation 



with the regulator held fixed. 

At one loop, the matter terms in the IPI effective action are 



(4.49) 



where 



(4.50) 



(4.51) 



Here, A is the ultraviolet cutoff. Invariance under the transformation Eq. ( 4.49 ) allows 
us to conclude that Fipi depends on 5*0 + 5*0 only in the invariant combination 



'S'o + >S'o — -— ^ In Zrfi. 



(4.52) 



This allows us to infer the 2-loop dependence of the matter kinetic term in Fipj from 
Eq. ( [4.51| ). We can then obtain 



327r 



4 \ ^ ^r'li'rfi 



(4.53) 



From this, we can read off the 2-loop RG equation for the soft masses arising from 
gauge interaction with other soft masses: 
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(4.54) 



(Note we have not specified a definition for the renormalized gauge coupling, but the 
result is invariant under changes of scheme for the gauge coupling.) 

If the gauge group contains a f/(l) factor, there is an additional contribution to the 
RG equation for the scalars from an induced Fayet-Illiopoulos term. In superspace, a 
Fayet-Illiopoulos term can be viewed as a "kinetic mixing" between the U{1) gauge 
field and that of the anomalous U{1) symmetries for the various matter fields. Note 
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that in the presence of bare soft masses, there is no symmetry forbidding such a term, 
so we have an addition contribution to the bare lagrangian 

5Cq = j<fe\Kr^QWiWrfl + h.C., (4.55) 

where Wi is the U{1) gauge field strength and 

(>V.,o)a = -^D^DJnZr,o = e^ml^ (4.56) 

is the field strength of the anomalous U{1). Eq. (|4.55 ) contains a linear term in the 
f/(l) auxiliary gauge field Z)i, forcing (-Di) 7^ and giving an additional contribution 
to the scalar mass. It is the running of this contribution that we now compute. 

The Fayet-Illiopoulos term is renormalized at one loop, and we obtain 

^^ 

Combining this result with the 1-loop renormalization of the matter wavefunction 
given in Eq. ( ^.51| ), we obtain an induced vacuum expectation value 



ipi 



j (fe I Lr,o + In j WiWrfl + h.c. + finite. (4.57) 



{D) = [y. <lr<o + E ^.CUrml,)^ In ^, (4.58) 

where the sum on J runs over the factors of the gauge group, and g,. is the U{1) 
charge of the field r. From this we can read off an additional contribution to the RG 
equation for the soft mass: 
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Recall that in Sect. 3 we showed that the RG equations for the soft masses above 
correspond to DR'. The present derivation shows that these RG equations follow 
as long as the theory is regulated and subtracted in a supersymmetric fashion. To 
further amplify this point, we give an illustrative application of these methods where 
we compute a soft mass as a finite calculable effect. 

Consider a toy model with bare lagrangian 



£0= jSeSaii{WWa)+\^.c. 



(4.60) 
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where r = 1, 2 are two copies of the same gauge representation. Suppose that the 
messengers (5i,2 have bare soft masses given by 

Z2,o = Zg,o [l + eWml] . (4.61) 

With this choice, the full theory has Str = 0, where Ai is the full mass matrix 
of the fields in the theory. However, if Mi ^ M2, the effective theory below the scale 
Ml has nonvanishing mass supertrace. The value of this supertrace is therefore a 
calculable effect in this theory. 

We could use the RG equations derived above to compute the soft masses in the 
low energy theory. We present here an alternative derivation of the supertrace that 
clarifies the methods used above. We assume M2 <^ Mi, and compute the q soft mass 
in the low-energy theory below the scale M2. With the choice of parameters made 
above, we can write 

Zi,o = e^°, ^2,0 = e-^°, Uo = -OWml (4.62) 

We can view Uq a,s a, "gauge" field for a single U{1) under which Qi and Qi have 
charge +1, Q2 and Q2 have charge —1, Mi has charge —2, and M2 has charge +2. 
Moreover, this U{1) symmetry is anomaly free, so we do not have to appeal directly 
to a Wilsonian picture of the anomaly. 

We now integrate out Q and construct the effective lagrangian below the scale 
M2. This has the form 

C" = y^^^Z^Ve^'" 5 + gauge terms, (4.63) 
where the U{1) symmetry enforces 



Z'' = /(|Mi|e-^MM2|e^°). 



(4.64) 



We can determine the function / by 1-loop running and tree-level matching to a scale 
At < M2: 
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where b {b') [b"] are the beta function coefficients in the full theory (effective theory 
below Ml) [effective theory below M2]. Using the 1-loop expressions for the gauge 



coupling g, and making the substitution |M 
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The first term corresponds precisely to the running of the soft mass between the 
scales Ml and M2, and the second term to the running between M2 and /i. There 
is no contribution from above the scale Mi because the contributions from the two 
messengers cancel. 



4.6 "Mediator" Models 



We now consider GMSB models where SUSY breaking is communicated less directly 
to the observable sector. We find that very generally in such models, the gaugino 
screening mechanism described in Sect. implies the gaugino mass is suppressed 
compared to the scalar masses by more loop factors than suggested by a naive analysis. 

We consider the "mediator" models introduced in Ref. [|]. We suppose that a 
SUSY-breaking sector communicates SUSY breaking to vectorlike fields Q and Q. 
The fields Q and Q are not charged under the standard-model gauge group. Rather, 
they are in a vector-like representation of a "mediator" gauge group Gmed- The 
connection to the observable sector is made through a vector like pair of fields T and 
T that are charged under both the standard- model gauge group and Gmed- These 
fields have a supersymmetric mass term My in the lagrangian, which may be the 
result of a dynamical mechanism H]. The lagrangian of this theory is 



c" = d^e 



(Q) - 



(4.67) 
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Af^Tf + 5L,tr«,,) + 5^Vtr(W^SM) + h-c 



+ (5£(g,Q,...), 

where 5£ contains the interactions that break SUSY. 

The holomorphic standard-model gauge coupling below the messenger scales M 
and the scale M^ is given exactly by 



5'sM(/i) - 5'gM(/io) + 
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(4.68) 



where h'^yi and 6sm are the standard-model beta function coefficients in the effective 
theory with and without the field T, respectively. This is independent of M, so 
the leading contribution to the gaugino mass comes from the In^^^ term in the real 
effective gauge coupling -Rsm, see Eq. (|4.19|) . The leading M-dependent contribution 
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to Zr arises at 4 loops, so the gaugino mass arises at 5 loops in this model, as opposed 
to the estimate of ref . 0] . Since scalar mass-squared terms arise at 4 loops, the gaugino 
mass is suppressed compared to the scalar masses in this model, posing a fine-tuning 
problem. 

To make this argument concrete, and to illustrate the power of our techniques, 
we explicitly compute the gaugino mass in the case where SUSY breaking is com- 
municated to the fields q and q by the vacuum expectation value of a singlet field 
X: 

(5£ = Jd^e XXQQ + h.c, (4.69) 

with (X), (Fx) 7^ 0. The reader uninterested in details can skip the remainder of this 
Subsection. 

We will do the calculation for the case where 

M = X{X) > Mt. (4.70) 

We further assume that Gmed is weakly coupled and unbroken down to the scale Mt- 
Below the scale M, the light fields are T, X, Qr, Vmed, and Vsm, and the effective 
lagrangian C consists of the terms in Eq. ( [4. 67] ) that depend on these fields. Below 
the scale Mt, the only light fields are X, Q^, Kned, and Vsm, and we denote the 
effective lagrangian by C 

Both the scalar and gaugino masses can be read off from Zr, the wavefunction 
renormalization factor in the low-energy theory. We therefore compute 



\nZ,{^)= r ^4ir{^')+ C %ir{^')+ [ ^7.(/^'), (4.71) 



where il'r) [1"] denotes the anomalous dimension in the theory C (C) [£"]; and 
M (Mt) is the matching scale at the mass of q (T), defined similarly to Eq. (|4.21 
For example. 
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where we have displayed the dependence on Z^ required by the "anomalous U{1)" 
invariance. This is important, because Z!p depends on M at 2 loops, giving the leading 
M dependence of the anomalous dimensions. We have 
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where 

ain|X| (87r2)2 {S'^^ + S'iy,fdhi\X\- ^ " > 

{T is not a light field in £, so there is no contribution from scales below Mt-) We 
therefore have 

d\n\x\ -{sn^yJM Pi' isu^^') + sU^^')r d\n\x\ ■ ^ - > 

(We see that the M-dependent part of Zr is independent of the renormalization scale 
/I.) The dependence of Z!p on the messenger threshold is identical to the calculation 
in GMSB, and we obtain 



(?ln|X| (Stt^)^ Jmt H Jh II' 

Prom this, we can obtain the gaugino mass 

/ N ^sm(/^) {Fx) ^ Tr d\nZr{ii) 
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Notice that the result scales like rux = ag^al^^^hi^ M/Mt, indicating that two loops 
are accounted for by 1-loop evolution. 



5 Effects from Other Thresholds 



Up to now, we have been focusing on effects that can be computed from the depen- 
dence on the messenger threshold. However, there are interesting models with other 
thresholds than can give rise to important SUSY-breaking effects in the low-energy 
theory. In this Section we analyze some illustrative examples. 



5.1 Flat Direction Effective Potential 

In the limit where SUSY is unbroken, the minimal supersymmetric standard model 
has a large space of fiat directions, directions in field space where the classical potential 
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vanishes identically. (For an exhaustive list, see Ref. @].) All of these flat directions 
will be lifted by SUSY breaking, and we are interested in computing the effective 
potential far out along one of these flat directions. For GMSB, the effective potential 
can be evaluated from 2- loop component diagrams such as those evaluated in Ref. [pT| , 
with the motivation of studying the cosmology of these fiat directions. We will show 
how to compute the effective potential without evaluating loop diagrams. 

We will explain our technique using a toy theory with an "observable sector" 
consisting of a U{1) gauge theory with Ng pairs of chiral fields q and q with charges 
+ 1 and —1, respectively. These are coupled to a "messenger sector" consisting of Nq 
pairs of chiral fields Q and Q and a singlet field X that parameterizes SUSY breaking. 
The lagrangian is 



C" = jd^e \z'l(q^e^q + q^e-^q 



+ Z'' (Q^e^Q + Q^e~^Q) + Z^XtX 



(5.1) 



+ jd'd^S"W''Wa,+\l.C. 

+ jd^9\Xqq + h.c. 

Even though X is a background field, we must include a "kinetic" term for X to 
account for the anomalous dimension of operators that depend on X. (This operator 
is just the contribution to the cosmological constant.) 

This theory has a single classical fiat direction with (g) = (g) . We want to compute 
the effective potential for (g) = (g) ^ {X). In this case, the largest threshold in the 
theory is at the scale 



M^ = g{M,)\{q)\, 



(5.2) 



where g is the U{1) gauge coupling. At this scale, the U{1) gauge group is completely 
broken. The fields that are light below this scale are Q, Q, and the fiat direction 
q = q, parameterized by a field Y defined as 



q={q) + Y, q= (g) + Y. 



(5.3) 



The background field X is also present in the low-energy theory. The effective la- 
grangian below the scale Mi is therefore 



£' = jd^e [z'q (Q^e^Q + Q^e-^Q) + Z'^X^X + Z'yY^Y 
+ Jd^9XQQ + h.c. + ---, 



(5.4) 
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where the ehpses denote higher-dimension operators. 

The next threshold of interest is the messenger threshold at M = A(X). Below 
this scale, the effective lagrangian contains only the fields X and y, and it is given 

by 



ZxX^X + ZyY^y] + • • • . (5.5) 



We are interested in the effective potential for Y in this effective lagrangian. When 
we continue the couplings into superspace, there will be contributions to the effective 
potential for Y from the Y dependence of Zx as well as the X dependence of Zy- The 
field Y docs not have renormalizable interactions below the scale Mi, so Zy does not 
depend on X at the renormalizable level. The contribution to the effective potential 
we are interested in is therefore 

V^^{\Y\)^-\{Fx)?Zx{\Y\). (5.6) 

We compute Zx using tree-level matching and 1-loop running. Using the RG 
equations 



dlnZ^ _ Nq dlnZ'x _ Nq y 

^ df, ^n^Z'J^Z'^' ^ df, ^n^Z'xZ^' ^'''^ 

we obtain 



_ ^/// NqX^ f^^ dfi 1 NqX^ r^' d^ 1 

Zx-Zx{iio)-^ jz^)~^^]M,Jz^y ^^-^^ 

where /io is a fixed renormalization scale used to define the theory. Note that Zx is 
independent of renormalization scale. Since we are interested in the Y dependence, 
we compute 

dZx _ NqX^ fMdfi 1 dlnZ'Qii^) 



d\n\Y\ 47r2 Jm^ fi Z'^{fi) dln\Y\ ' ^^'^^ 
Zq does not run in the effective theory C, so we have Zq{ii) — Zq{Mi), which gives 

dlnZ'^ifx) _g\M^) 



d\n\Y\ 87r2 

In this way, we obtain 



(5.10) 



' (47r2)2 Z2(Mi) M' ^ ' 

Note that Mi depends on so this result automatically gives the RG- improved 
form of the effective potential. 
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5.2 (S)axion Potential 



There are a number of models for physics beyond the standard model that involve 
the spontaneous breaking of a global symmetry at large energy scales. For example, 
"invisible" axion models invoke the breaking of a global ?7(l)pQ symmetry at scales 
10^°-10^^ GeV in order to solve the strong CP problem. Other global symmetries 
that may be spontaneously broken include lepton number and flavor symmetries. 

The breaking of a global symmetry will give rise to a massless Nambu-Goldstone 
boson (NGB) for every broken generator. If the global symmetry is broken at a 
scale where SUSY is (approximately) unbroken in the visible sector, then the light 
bosons must form complete chiral supermultiplets. There are therefore extra scalars 
whose mass is protected by SUSY.0 We call these fields SNGB's. The SNGB fields 
parameterize non-compact directions in the vacuum manifold in the limit where SUSY 
is exact, and different points along the flat direction correpond to different values for 
the scale at which the global symmetry is broken. The SNGB fields will acquire 
a potential after SUSY breaking, which determines the vacuum expectation values 
along the fiat direction. 

As an example, we consider an axion model with colored fields R and R whose 
mass is determined by the vacuum expectation value of a field $. If we write 

$ = ($) + A, (5.12) 

the imaginary part of A is the axion, while the real part is the SNGB. The lagrangian 



IS 



+ Jd^9K(l>RR + h.c. + ---, 



(5.13) 



where we have omitted the messenger sector and standard-model fields, see Eq. (^.1|). 
The fields R, R therefore have a mass 

Mr = ^P--. (5.14) 

Below this scale, the effective lagrangian C is simply that of ordinary GMSB together 
with a kinetic term for the field $ (see Eq. ( [4.1|) ). Below the messenger threshold M 
the effective lagrangian C is that of the standard model together with kinetic terms 

^If a non-abelian symmetry is broken, some of the Nambu-Goldstone bosons can belong to the 
same chiral supermultiplet, but it can be shown that there are always some "extra" scalars. 
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for the singlets $ and X. The wavefunction parameter Zx in this effective lagrangian 
depends on the R mass, and this contains the leading contribution to the effective 
potential for the saxion field. 

We can compute Zx using 1-loop running and tree-level matching: 

Zx = Z'^if^o) + ^ + r - 2'xif^hM, (5.15) 

JflO fl JMr fl 

where 7x is the anomalous dimension of X as defined in Eq. ( ^.7|) . The parameter 
Zx does not run in this effective theory, so we need not specify a renormalization 
scale for it. We compute 

dZx r*^ dfi d 



(5.16) 



ain|$| 7m, mn|$|[^^^^)^^(^)] 

Nq\\\^ fM dii 1 51nZ'(/i) 



47r2 J Ma /i Z'^{ii) (91n|<l>| 
The right-hand side is evaluated using 



91n|<l>| J Mr /i' 91n|$| 



Cq f^' rf// ,4 , d 
,9 yf' )' 



Atc^Jmr /i' '"^'Slnl^l \g'^{fi') 



which gives 



U 2.2 r ^9"if^^ (5-17) 
(47r^)^ JMr /i' 



As before, this gives the RG-improved form for the effective potential. Note that 
the slope of the potential is negative, indicating that the saxion vacuum expectation 
value is driven away from the origin. 

In the opposite limit Mr -C M, it is easy to see that the potential also decreases 
as a function of Mr. In the effective theory at the scale M, R and R get a positive 
soft mass-squared from GMSB while $ has zero soft mass. However, the Yukawa 
coupling K^RR drives the $ soft mass^ negative in running between M and the 
scale Mr, where R and R are integrated out. (This contribution is analogous to the 
negative contribution to the Higgs mass-squared from the top Yukawa coupling.) 

Thus in all regions, the potential prefers to push the saxion vacuum expecta- 
tion value, and hence the axion decay constant, to larger values. Therefore new 
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interactions are needed between the axion and GMSB sectors in order to stabilize 
the axion decay constant in the cosmological and astrophysically desirable window 
between lO^O-lQi^ GeV. 

6 Conclusions 

In this paper, we have shown that the renormalization of soft SUSY-breaking terms is 
completely determined by the renormalization of SUSY-preserving terms if the regu- 
lator is supersymmetric. This allows us to calculate certain SUSY-breaking effects in 
gauge-mediated theories by performing a supersymmetric calculation and "analyti- 
cally continuing" the result into superspace. The method is very powerful, and allows 
the calculation of interesting effects at 3-loop order and higher by purely algebraic 
manipulations. 

The formal results that justify these calculations are easy to state in superspace 
if the soft SUSY breaking terms are parameterized by ^-dependent terms in the su- 
persymmetric couplings. If the theory is regulated in a supersymmetric manner, then 
SUSY is formally preserved if we regard the bare couplings as superfield spurious. 
Our result is that there is a definition of the renormalized couplings that can be 
similarly grouped into supermultiplets. Specifically, the renormalized couplings Kji 
are related to the bare couplings Kq via a superfield relation of the form 

KM^f{K,,A,f,). (6.1) 

The function / determines the renormalization of the supersymmetric couplings as 
well as the soft SUSY breaking terms, and is the basis for the analytic continuation 
into superspace. An analogous relation holds between the (renormalized) couplings 
of an effective theory and the couplings in a more fundamental theory. 

This leads naturally to a definition of the renormalized gauge coupling chiral 
superfield 

SW = ^-^-»=^ (6.2) 

as a holomorphic object that is renormalized only at one loop (to all orders in per- 
turbation theory). However, the subtraction that defines S{fj.) is not invariant under 
constant rescaling of the fields, so the components of S{fj.) do not correspond directly 
the usual renormalized couplings. The real superfield 



R^S + S^-^ HS + S^)-Y^^\nZr + 0{{S + S^)-'). (6.3) 
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is invariant under the field rescaling. (Here, r runs over the matter representations 
of the gauge group G and Tj. is the index of r; is the wavefunction factor for the 
fields in the representation r.) We show that the lowest components of R 

^ -R{f^)\, ^^=m\o^, (6.4) 



are precisely the IPI gauge coupling and gaugino mass defined by Euclidean sub- 
traction or by minimal subtraction in dimensional reduction. The 0{{S + 5'''')"^) 
corrections account for possible scheme dependence in the definition of R. Eq. ( |6.3| ) 
and much of the story leading up to it is very similar to the results of Refs. 0, but 
we emphasize that all quantities are finite renormalized quantities, and no reference 
is made to the Wilsonian renormalization group. 

The 6'^ 6"^ component R is given at lowest order by 



R\ 



(6.5) 



and governs the RG evolution of dimension-2 soft terms. In dimensional reduction 
R\0202 corresponds to a 1/e counterterm for the e-scalar mass. -R| 512^2 can also be given 
a IPI interpretation: it corresponds to a non-local correction to the propagator 
of the gauge supermultiplet. In the context of dimensional reduction and (modified) 
minimal subtraction, our results imply that the simple extension l/g'^{n) —>■ R{n) 
automatically picks out the so-called DR' scheme. 

In practice, this result allows one to simply compute the SUSY breaking compo- 
nents of R (for example) by computing the lowest component as a function of the 
supersymmetric bare couplings (or couplings in an underlying renormalized theory). 
This is a supersymmetric calculation, but taking ^-dependent components of the re- 
sult determines the low-energy SUSY breaking parameters. For instance, we have 
shown that the 2-loop RG equations for soft terms in DR' are directly derived from 
the supersymmetric /5-functions and anomalous dimensions. 

More remarkably, this approach can be used to relate leading-log effects computed 
using the renormalization group to finite effects, since the result of taking higher 6 
components of a logarithm gives effects that are not logarithmically enhanced: 



^ InM 



167r2 



1 M\e-e- 
0202 IGvr^ M\ 



(6.6) 



In this way, we can obtain finite SUSY breaking effects at high loop order from simple 
algebraic calculations. Models with low-energy supersymmetry breaking mediated by 
perturbative interactions are the natural arena to apply our method. Indeed, it is 
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precisely in these theories that it makes more sense to worry also about subleading 
RG evolution: this is because the boundary conditions for soft terms are in principle 
calculable with comparable accuracy. 

Our technique was used to compute a variety of effects at 2-loop order and beyond. 
We computed for the first time the complete subleading corrections to the gaugino 
masses (2-loop) and scalar masses (3-loop) in gauge-mediated models; we showed how 
to compute the effective potential for SUSY flat directions lifted by gauge mediation 
(2- and 3-loop). We also proved that gaugino masses are screened from higher- loop 
corrections involving couplings in the messenger sector. Therefore, in the standard 
gauge mediated scenario, gaugino masses are rather insensitive on details of the model. 
Moreover, this result also shows that if the gaugino masses are not generated at one 
loop (as in the standard case) they will be generated only from the light matter fields, 
and will generally be too light. This shows that gauge mediation is the unique way 
to generate scalar and gaugino masses of the same order through loop effects. 
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